Flapping motion and force generation in a viscoelastic fluid 
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In a variety of biological situations, swimming cells have to move through complex fluids. Sim- 
' ilarly, mucociliary clearance involves the transport of polymeric fluids by beating cilia. Here, we 

consider the extent to which complex fluids could be exploited for force generation on small scales. 
We consider a prototypical reciprocal motion (i.e. identical under time-reversal symmetry): The 
Q , periodic flapping of a tethered semi-infinite plane. In the Newtonian limit, such motion cannot be 

, used for force generation according to Purcell's scallop theorem. In a polymeric fluid (Oldroyd-B, 

I and its generalization), we show that this is not the case and calculate explicitly the forces on the 

' flapper for small-amplitude sinusoidal motion. Three setups are considered: a flapper near a wall, 

r"j^ ' a flapper in a wedge, and a two-dimensional scallop-like flapper. In all cases, we show that at 

O . quadratic order in the oscillation amplitude, the tethered flapping motion induces net forces but no 

average flow. Our results demonstrate therefore that the scallop theorem is not valid in polymeric 
fluids. The reciprocal component of the movement of biological appendages such as cilia can thus 
^ ' generate nontrivial forces in polymeric fluid such as mucus, and normal-stress differences can be 

' isJ ' exploited as a pure viscoelastic force generation and propulsion method. 

; 

^ ' PACS numbers: 47.57.-s,47.15.G-,47.63.-b,47.63.Gd 

^ ! I. INTRODUCTION 

i> : 

^ . 

, Locomotion is a subject with a rich history, and with relevance to all length scales in biology, not only in the animal 
i kingdom - from large cetaceans to spermatozoa - but also to bacteria, protozoa, and algae l], Q . 
I A particularly active area of research concerns cell motility in viscous fluids. Classical work in the 1950-70's have 
00 , elucidated the fundamental principles of fluid-based propulsion on small length scales , but a number of 

physical issues are still unresolved. Most notably, the extent to which its surrounding environment plays an important 
role in the swimming of an individual cells is not well understood. 
^ . One relevant example is the locomotion of cells in complex fluids, such as the swimming of spermatozoa in cervical 
5^ , mucus 0, Q| . Cervical mucus is a high- viscosity cross-linked polymeric gel [lo| located at the entrance of the uterus 
along the female reproductive tract. Its most notable feature is to possess relaxation time scales in the 1-100 seconds 
range which are much longer than the typical time scales involved in the beating of flagella 



- '—I 

X 



(typical frequencies in the tens of hertz). As a consequence, the mechanics of locomotion in cervical mucus is 
significantly different from locomotion in Newtonian flows. Experimentally, spermatozoa swimming in complex fluid 
display a change in their beat pattern (smaller amplitude, wavelength), an increase in their beat frequency, and (as 



a result) a change of their swimming kinematics IT], [iSl, |l9|, |20|, |2l|, [22, l23| . Theoretically, recent work predicted that 
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cells actuating their flagella in a wave-like fashion should systematically swim slower in a polymeric fluid than in a 
Newtonian fluid 
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From a general standpoint, however, there is no reason to expect that going from a Newtonian to a viscoelastic 
fluid should systematically lead to a degradation of the swimming performance. Since a viscoelastic fluid displays 
nonlinear rheological behavior [27|, |28|, |29|, |30| , in general it should be possible for cells to exploit these nonlinearities 



to gain propulsive advantage. The goal of this paper is consider the simplest possible body deformation able to take 
advantage of these nonlinearities for force generation. 

In a Newtonian flow, it is known that there exists a class of body deformations - those where the sequence of body 
shapes are identical under time-reversal symmetry, called "reciprocal" - that lead to exactly zero net propulsion in the 
Newtonian limit in the absence of inertia 
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Such result is known as Purcell's scallop theorem [31l |. 
and is due to the linearity and time-reversibility of Stokes equation of motion. The prototypical reciprocal motion 
has only one degree of freedom, such as a flapping back-and-forth motion (e.g. the flapping of a wing, or the motion 
of a scallop). As a result of the constraints of the scallop theorem, swimming microorganisms are observed to always 
deform their body or their flagella in a wave-like fashion in order to move. Similarly, mucociliary clearance in our 
lungs involve fluid transport by arrays of cilia beating periodically in a non-reciprocal manner [3^ . 

In this paper, we consider tethered flapping motion in a variety of two-dimensional settings in the absence of 
inertia as the simplest model of biological reciprocal motion. We calculate the flow fleld perturbatively in the flapping 
amplitude for a polymeric fluid (Oldroyd-B and generalized models). Although no net flow is obtained at second 
order for a periodic actuation, we show that net forces are generated. Reciprocal motions arc therefore able to 
harness normal stress differences to generate net forces, and the scallop theorem breaks down in polymeric fluids. An 
implication of these results is the possibility for the reciprocal component of biological appendages, such as cilia, to 
generate nontrivial forces when moving in polymeric fluids such as mucus. 

In many ways, our calculations are reminiscent of classical work by GI Taylor and HK Moffatt on similarity solutions 
for the Newtonian paint scraper problem 
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39( 1 . and two-dimensional viscous vortices 



4C , 41 1 . They are also 



related to the study of polymeric stress on steady (or "acoustic" ) streaming [42l . |43| . l44l |45| . |46| . |47| . |48| . |49| . ISOl . l5ll | , 
which is the physical phenomenon by which the periodic actuation of a solid object is able to exploit inertial forces to 
generate net flows. In particular, polymeric stresses lead to a reversal of the net flow direction 42, 143], |4J, |45( . Here, 
we are concerned however with purely viscoelastic forces, and ignore the influence of inertia. For work related to the 
current study in the limit of weakly viscoelastic flows (small Deborah number), see Ref. [5^. 

The paper is organized as follows. The general setup for our two-dimensional flapping calculations is presented in 
^Jn] We then calculate in detail the time-averaged forces on a semi-infinite flapper near a wall, and more generally 
in a wedge, for an Oldroyd-B fluid in mill The forces on a two-dimensional scallop are next calculated in mVl We 
generalize our results, and discuss their implications for biological propulsion in ^ 



II. VISCOELASTIC FORCES IN SMALL-AMPLITUDE FLAPPING: GENERAL RESULTS 

A. Setup 

We consider a semi-infinite two-dimensional plane flapping with small amplitude in a viscoelastic fluid. The three 
setups we will consider are illustrated in Fig. [T] (a) a flapper near a flat wall, (b) a flapper in a wedge, and (c) a 
scallop-like body made of two symmetric flappers. We will first calculate the flow on one half of the flow domain (e.g. 
< 6* < a in Fig. [1^) and deduce the flow in the second part of the domain (e.g. a < 6 < n in Fig. [1^) by analogy. 
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FIG. 1; (color online) Cases considered in this paper: (a) two dimensional flapping motion near a no-slip surface (flapping of 
amplitude e around a mean angle q); (b) flapping motion in a no-slip wedge (wedge angle, 2/3); (c) scallop- like flapping motion 
(mean angle 2q between the two arms of the scallop). In all cases, the hinge point of the flappers are flxed in the lab frame. 

B. Constitutive equation 

We assume in this paper incompressible flow in the absence of inertia. The velocity field is denoted u and the 
pressure field p. In the fluid domain, Cauchy's equation of motion is therefore written as 

V-u = 0, (la) 
Vp = V-T-, (lb) 

where t is the deviatoric stress tensor. We flrst consider polymeric fluids described by the Oldroyd-B constitutive 
relationship 27, 28, 3 



Ai x= ?7[7 + A2 7], (2) 



with 7 = Vu + Vu"^. In Eq. ^ we have defined, for a tensor A, the upper-convected derivative A= ^ + u • VA — 
(Vu-'" • A + A ■ Vu); Ai and A2 < Ai are the relaxation and retardation times of the fluid respectively. We consider 
periodic flapping motion with frequency uj. We non-dimensionalize rates and stresses in Eqs. ([IJ-iH]) by co and yyw 
respectively, and lengths by some (arbitrary) length scale along the flapper. The dimensionless equations are therefore 
given by 

V-u = 0, (3a) 
Vp = V-T, (3b) 

V ^ 

T + Dci T = 7 + De2 7, (3c) 

where Dei = Aiw and De2 = A2ti-' are the two Deborah numbers for the flow, and where we use the same symbols 
for convenience. The Oldroyd-B constitutive relationship, which can be derived analytically from a dilute solution of 
perfectly elastic dumbbells [30], correctly models polymeric fluids up to order one Deborah numbers. For larger values 
of De, more accurate models are necessary to capture the correct extensional rheology of polymeric flows. One of such 



models is considered in the discussion (FENE-P), and we obtain the same results (see also 2J|). Other Oldroyd-like 
models with more complex rheological characteristics are also discussed in fJV] 
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FIG. 2: (color online) Setup and notation for the flapping calculation: the inflnite flapper (black) has an opening angle which 
oscillates in time around a mean value, a (red). For a point in the fluid domain at a distance r from the hinge point, the polar 
vectors are denoted er{9) and eg{9) (blue). The polar vectors attached to the average flapper position are denoted (t) and (n) 
(red). 

C. Streamfunction 

Given the geometrical setup, it is most convenient to use polar coordinates, and introduce the streamfunction, 
defined by Ur — {d'i'/d9)/r, ug = —d'^/dr. The polar vectors are denoted e,.(6') and 6^(6*) (see Fig. [5]), and the 
velocity field is denoted u = Ur&r + ugeg. 

D. Perturbation expansion 

The calculations in the paper are done perturbatively in the amplitude of the flapping. We denote the instantaneous 
position of the flapper as 9 — a+eQ(t), with &{t) = cost an order-one oscillatory function, and perform the calculation 
in the asymptotic limit e ^ 1 using domain perturbation. We look therefore to solve the problem using an expansion 
of the form 

{u,5',T,p,cr} = e{ui,^'i,Ti,pi,cri} (4) 

+ e^{u2,*2,T2,P2,Cr2} + 

where cr = —pi + x is the total stress tensor, and where all variables in in Eq. (j4]) are defined in the averaged domain 

o<e <a. 

E. Instantaneous vs. average positions 

Since we perform a domain-perturbation expansion, we have to pay close attention to the distinction between 
instantaneous and average geometry. We use planar polar coordinates. The vectors and eg on the flapper are 
functions of the azimuthal angle which oscillates in time: er{a + eQ{t)), eg{a + eQ{t)). In the paper, we use the 
notations (t) = 6^(0;) and (n) = eg{a) to denote the average directions along and perpendicular to the flapper axis. 
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F. Fourier notations 

It is easier to solve for the flow field using Fourier notations. The actuation is defined as 8(t) = 5R{e'*}, Q{t) = 
5R{ie**}. Since quadratic nonlinearities arise from boundary conditions and the constitutive model (see below), we 
anticipate a Fourier decomposition of the velocity field of the form 

ui = 5R{uie**}, (5a) 
U2 = 5R{iif +a^'^e2^*}, (5b) 

and we use similar notations for all other vector and scalar fields in the problem. 



G. General equations 



We now look in detail at the first two orders of Eq. ([3]). At leading order, the nonlinear terms disappear, and we 
obtain 

rr+De,— ^7,+De.^, (6) 



which becomes, in Fourier notations, 



1 + iDe2 ~ 
f + zDci 



We then take the divergence of Eq. ([Tj) and the curl of Eq. (|3bp . to obtain the equation for the as streamfunction 

d_ 

dt 



1 + I-DC2 ) V**i = 0. (8) 



At order 0(6^), Eq. ^ becomes 



1 + ^Dei) T2-(l + I-DC2 ) 72 (9) 



dt 'J ^ \ dt 
= De2 [ui • Vii - (*Vui • 7i + 7i • Vui)] 

-Dei [ui • V-Ti - (*Vui -Ti+Ti ■ Vui)] 
Using Eq. the average of Eq. ([9]) becomes, using Fourier notations, 



-^(0) _ - 1 De^^Dei . ~ 



-(*vii^ -^1 + ^1 - vat)] , 

which can be evaluated using knowledge of the solution at order 0{e). 



H. Boundary conditions 

Associated with the equations derived in §11 Gl and §11 1[ we have to obtain boundary conditions at each order. On 
steady surfaces, the boundary conditions can be simply stated. For the steady no-slip surfaces in Fig. [1^ and Fig. [1)3 
we have u = on the boundary. For the symmetry plane 9 = in Fig. [TJ: we have the no-penetration conditions 
U0 = 0, as well as symmetry, dur/d6 — 0. 
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The other boundary condition occuring on the flapper is the statement of its rotation, u(r, 9 = a + e8(<)) = rVl{t)eg 
(with ri[t) ~ 60), which is obtained order by order by performing a domain perturbation expansion. Consider a general 
vector field h(r, 6*). A Taylor expansion for h on the flapper position gives 



h(r, 6* = a + eG) = h(r, a) + eO _ 

Ou 



+ ..., (11) 

(r,a) 



which we apply for h = u — reQee = 0, and therefore obtain, at order by order, 

ui = re(n), (12a) 

U2 = -e^-ree(t), (12b) 

when evaluated ai — a. In Fourier notations, we use Eq. ([5]) into Eq. (|12p and obtain the boundary conditions for 
the first two Fourier components as given by 

Ui — ir{n), (13a) 
- ^2W 



evaluated at 6 = a. As we see below, we do not need to solve for the Fourier component u. 



r.(2) 



I. Time-averaged stress on the flapper 

We are interested in calculating the time-averaged stress on the moving flapper. The unit normal into the fluid is 
—eg and therefore the instantaneous stress is given by 

5f = -(Treer - 0-0060. (14) 

Using Taylor expansions for both the stress components and the direction vectors in Eq. Ijl4p (i.e. applying Eq. (jlip 
for h = —cr • ee), we find that there is no average stress at first order on the flapping surface, while the stress at order 
O(e^) along the averaged directions of the flapper, i.e. ((t), (n)), is given by 

(t) • 5f2 = -<y2.re - eit) - ai,,^) , (15a) 

(n) • 5h = 'O^M - m + ^i.r^ , (15b) 

evaluated at 9 = a. In Fourier notations, and given that 6(t) = cost, they become 

(t) • {5h) = -5R + \ - o,,eo^ } , (16a) 

(n) • (^f2) = [of}, + \ + a.^re^ | , (16b) 

evaluated at 9 — a. Thus, in order to compute the average stress at this order, we only need to know the Fourier 
components of the flow Ui and iif, and do not need to solve for the time- varying contribution given by U2^''. 

Finally, on the no slip walls, there is no moving boundary. The average stress is therefore simply given there by 
the Fourier component, a-f, evaluated on the wall. 
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III. SINGLE FLAPPER 

We now solve for the case where there is only one flapper. As stated above, we solve the flow in the domain 
< 9 < a, and deduce the total stress by analogy. 

A. Solution at 0{e) 
1. Velocity field 

Taking the Fourier transform of Eq. ([5]), we get V^^i = 0. With the boundary conditions given by Eq. (|13a[) . we 
obtain the solution at first order 

~ ir^(— sin26' + tanacos2^ + — tana) 

^ " 4(tana-a) ' ^ ' 



which is the same as for the Stokes problem 



0. 



2. Pressure 



The deviatotic stress is given by Eq. ([7]), so we have 

1 + iDei J (tan a ~ a 



- I \^ 1 .T, (18) 



with 



r(l,l) = 1 - cos 26* - tan a sin 20, (19a) 

T(f,2) ^ sin 20 - tan a cos 26*, (19b) 

r(2,l) = T(l,2), (19c) 

r(2,2) = -r(l,f). (19d) 



The pressure is found by integrating Cauchy's equation, Eq. ([3b|) . We calculate the divergence of the deviatoric stress 
tensor, Eq. ((TS]), and integrate it to find 

Pi = L ^ -T, 7 - ' (20) 

\ i + iDci / tan a ~ a. V ^'o / 
which has an integrable singularity at the origin (ro is a dimensionless microscopic cutoff length) j4o| . 

3. Contribution to second-order stress 

With the stress calculations above, we get at the average flapper position 

d&i^re -. 1 _ Dei - De2 / 2 

89 '^^■fis j 1 + De^ \ tan a — a 



^\-^~-^..ee\ = ^^^^^(rrr^) (21a) 



1 + In ( - 



daigg ^ 1 Dci — Dc2 / tana 
89 ^i^'-fj 1 + De^ \a-tana 

which, as expected, are equal to zero in the Newtonian limit, Dei = De2 
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B. Solution at 0(e^ 



1. Velocity field 



a. Boundary conditions. The boundary conditions for the average flow at the location of the flapper, Eq. (|13b[) . 
become 



.-.(0) 



u 



-(t). 



^ 2(tanQ; — a) 
h. Average flow. FoUowing Eq. (jlO[) . we calculate 

ui • V^i - (* Vui • ^1 + i'l • Vui) = 

with 



A 



(tan a — a) 



2 + tan^ a - 2 cos 26* - 2 tan a sin 20 
— (tan a — 26) (sin 29 ~ tan a cos 2ff) , 
-sin26'(l + tan^ a) 
+26'(cos 29 + tan a sin 26), 
A{2,1) = Ail,2), 

2 + tan'^ a — 2 cos 26 — 2 tan a sin 26 
+ (tan a - 26) (sin 26* - tan a cos 26*) . 



A(l,2) 



A{2,2) 



Since Ui is purely imaginary, the equation for the zeroth harmonics becomes then 



f (0) _ ^(0) 



1 De2 - Dei 

2 (tana -a)^ 1 + zDci 



In order to obtain the streamfunction at second order, we note that V • A = to get 



Combining with Cauchy's equation, Eq. (j3b|) and taking the curl, we get 

V4^^°) = 0. 



(0) 



(22) 

(23) 

(24a) 

(24b) 

(24c) 
(24d) 



(25) 



(26) 



(27) 



Given that the boundary condition, Eq. (j22p . is purely imaginary, we obtain therefore the following result: there is 
not average flow at order (^2*''' = 

2. Pressure and stress field 



We now have 



(28) 



and therefore the pressure at order is uniform (zero). 
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3. Contribution to stress 



/Dei — De2 



1 



On the flapper (9 = a), we have 



^2,00 V 1 + «Dei J tana - a' 
whereas on the bottom wafl {0 — 0) we obtain 



2'''^ V l + «Dei J a -tan a' 
(0) / Dei — De2 ^ tan a 



-(0) 



0) - 0, 
0) 0. 



(29a) 
(29b) 



(30a) 
(30b) 



C. Total time-averaged stress 

The calculations above allow us to evaluate Eq. p6p and obtain the time-averaged stress on the flapper at second 



order 



(t) • {Sh) 
(n) . {Sh) 



Dc2 — Dci / In {r/ro) 



1 + Dc^ \tan a ~ a 
De2 — Dei / tana 



2(1 + Def) Vtana 



(31a) 
(31b) 



D. Flapper near a wall 

We first apply our results to the case illustrated in Fig. [T^, that of a flapper near an infinite no-slip surface. We can 
assume a < 7r/2 without loss of generality. Combining Eq. (|3ip for angles a and n — a leads to the total time-averaged 
(dimensionless) stress on the flapper at second order 



7rDe(77,,/?7 - 1) 



tana(n) + 21n l;^- ) (t) 



2(1 -I- De^) (tana — a)(7r -I- tana — a) 



(32) 



Note that in Eq. ([32]) . we have replaced the two Deborah numbers by using the definition De = Dei, and ruis the 
solvent viscosity, rjs = 1^X2/ \i {ris < f]). Note also that the stress has an integrable singularity at the origin 4C|. 

The main conclusion from Eq. ([32]) is that the scallop theorem breaks down for a viscoelastic fluid: Although 
the flapper is acting on the fluid with a reciprocal (time-periodic) motion, it is able to generate net forces. Since 
'7s < Vi see that {6f2) is directed along the — n and — t directions. In the Newtonian limit, we have De = 0, and 
recover therefore the scallop theorem, {di2) = 0. Physically, the quadratic forces in Eq. ([32]) stem from normal-stress 



differences which arise in polymeric fluids due to flow-induced stretching of polymer molecules 
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281, 



29 
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E. Flapper in a corner 



The formula above for a single flapper can be generalized when the flapper is located in a wedge-like geometry of 
size 2/3 (Eq. [32] is the particular case /3 = 7i'/2). We can assume a < /3 without loss of generality. In that case the 
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total time-averaged stress along the flapper is given by 

De(77,/r;-l) 



(n) . {Si2) 



(t) ■ {Sh 



2(1 + De^) 
a tan(2/3 ~ a) + {a — 2/3) tan a 



(tana — a)(tan(2/3 — a) — 2/3 + a) 
De(7?,/77 - 1) 



1 + De^ 



In 



(33a) 



(33b) 



tana + tan(2/3 - a) - 2/3 



(tana - a)(tan(2/3 - a) -2/3 + a)_ 

The sign of the force in Eq. (|33p is easily obtained numerically. In particular, by tuning the values of /? and a, it is 
straightforward to show that all four sign combinations for the components (n) • ((5f2) and (t) • ((5f2) are possible. 



IV. A SCALLOP-LIKE FLAPPER 



We examine here the case illustrated in Fig. [TJ; for a tethered scallop-like two dimensional flapper. The difference 
with the calculations above lies in the boundary conditions, as explained in ^11 HI We perform the calculations below 
for the domain < 6* < a, and proceed to obtain the total time-averaged stress by analogy. 



~ (sin 261 - 261 cos 2a) 

*i = . „ . ■ (34) 



A. Solution at 0(e) 
1. Velocity field 

With the appropriate boundary conditions at 6 = 0, the solution at first order is given by 

ir2(sin26 - 26 cos 2a; 
2(2a cos 2a — sin 2a) 

2. Pressure and stress field 

The first-order deviatoric stress tensor is given by 

1 + iDi 

1 -I- iDei (2a cos 2a — sin 2a) 



1 + iDe2 2z _ , , 



with 



T(l,l) = cos 26 -cos 2a, (36a) 

r(l,2) = -sin 26, (36b) 

T(2,l) = r(l,2), (36c) 

T(2,2) = -T(l,l). (36d) 



The divergence of Ti can be integrated to lead to the pressure 

1 + iDe2 / 4?; 
1 + zDci I, tan 2a — 2a/ Vro 
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3. Contribution to second-order stress 



As part of Eq. (I16p , we have the foUowing 



^ f dd^.re . \ (Dci-De2)4[l + ln(r/ro)] 

\ de ] l + Del tan2a-2a ^ ' 

sa{ddifie , , 1 (Dei-De2) 2 tan 2a 

y 89 J 1 + Def 2a - tan 2a 



B. Solution at 0(e^ 



1. Velocity field 



a. Boundary conditions. Here, the boundary conditions from Eq. (|13bp become 



(0)^-2a + tan2a 

2 2 2a-tan2a^ ' ^ ' 



b. Average flow. After some algebra we get 



111 • Vi'i - ('Vdi • + ^1 • VGi) (40) 
4B 



(2a cos 2a — sin 2a)^ 



with 



(41) 



B{1,1) = 1 + cos 2a(26' sin 26* + cos 2a - 2 cos 26*) (42a) 

5(1,2) = cos2a(26'cos26'-sin26'), (42b) 

5(2,1) = 5(1,2), (42c) 

5(2,2) = 1 + cos 2q (cos 2a -2 cos 26* - 26* sin 26*). (42d) 

As in the case of a single flapper, we observe that V • B = to obtain that there is no average flow at second order, 
= 0. 

2. Pressure and stress field 

We have 



.(0) _ 2 (Dci - Dc2) B_ 

1 + iDei (2a cos 2a - sin 2a)^ 



as well as zero pressure {p2 = 0), so that the average stress tensor on the flapper is given by 

.(0) _ 2(Dei-De2) 1 

"^.re - i^^iDei 2a -tan 2a' ^ ' 



(0) ^ 2(Dei-De2) tan 2a _ 

2'^^ 1 + iDei tan 2a -2a ^ ' 
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C. Total time-averaged stress 

With the results of Eqs. (ISS)) and (|44|) . we obtain the time-averaged stress on the flapper as given by 



2(De2-Dci) ln(r/ro) 
1 + De^ tan 2q: - 2a ' 



/+\ /xf \ ^\i^^2 - ^J^i i^H'/'o; 

W-('5f2) = , , ^2 .„„o„. o„, ' (45a) 



I \ i<:r\ Dc2 - Dei tan 2a 

(n) • ((5f2) = — (45b) 

^ ' ^ ^' 1 + Dei tan2a-2a ^ ' 

D. Time-averaged stress on scallop 

Once we have the calculation for flapper in the domain < < a, we can use two reflection symmetries to obtain 
the time- aver aged stress on the entire scallop. We get 

47r cos a [2 In (r /ro ) — tan a tan 2a] 
(tan 2a - 2a) {2tt + tan 2a - 2a) ' 

Interestingly, in this simplified two-dimensional setup, the sign of the net stress can only be established when 7r/4 < 
a < 7r/2. Indeed, when < a < 7r/4, the numerator of Eq. ((46|) has an undetermined sign and depends on the size of 
the scallop. If we denote by ac the solution of tan2ac = 2(ac — tt) (ac w 51.27°), the force is in the —x direction for 
ckc < a < 7r/2 and in the -I- a; direction for 7r/4 < a < ac- 



V. DISCUSSION 



A. Viscoelastic propulsion. 



In this paper, we have used asymptotic calculations to derive the time-averaged forces acting on a tethered flapper in 
a polymeric (Oldroyd-B) fluid. The main results we obtained are Eqs. ((32|) . ([33| and ((46|) . which give the dimensionless 
stress acting on a flapper near a wall, in a wedge, and on a two-dimensional scallop respectively. 

The calculations presented above were made under a number of simplifying assumptions. First and foremost, we 
have assumed the flapper to be tethered, and therefore do not solve for the swimming motion that would be caused 
by the net forces it is generating. The second simplifying assumption is that of a two-dimensional geom etry. In that 



regard our work is the extension of classical work on similarity solutions for viscous corner flows 371, l38l |39| . The 
advantage of considering such an idealized geometrical setting is to allow us to solve rigorously the equations of motion 
for the viscoelastic flow, and obtain in a closed form the resulting time-averaged forces. 

Within the framework delimited by these assumptions, our results show rigorously that, using a purely harmonic 
motion, it is possible to escape the constraints of the scallop theorem in a viscoelastic fluid, even in the absence of 
inertia. Since the flapper is moving with a pure sinusoidal motion, the forces arising on the flappers do not originate 
from any rate-dependent mechanisms (which could be the case, say, if it was opening fast, and closing slowly). Instead, 
the net force arises from a rectification of the time-periodic actuation by normal-stress differences (i.e. flow-induced 
tension along streamlines due to the stretching of polymeric molecules). We note that since the flow is not viscometric, 
the explicit relationship between classically measured rheological parameters (such as normal stress coefficients) and 
the magnitude and direction of the time-averaged forces is not evident. A similar situation - but simpler to analyze 
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- arises in a cone-and-plate rheometer where the flow is viscometric. Under the oscillatory motion of the rheometer 
quadratic forces proportional to the first normal stress coefficient of the fluid act to separate the cone from the plate 
[2?! . These are the same forces which are responsible for the famous non-Newtonian rod-climbing effect . 



B. Order of magnitude and biological relevance. 

In our work, the hinge point of the flapper is fixed in place. The force arising from the fluid correspond therefore 
to tethered swimming, and are directly relevant to the motion of cilia-like biological appendages. In particular, our 
results show that, as a difference with a Newtonian fluid, the component in the motion of a cilium which is reciprocal 
(its periodic back and forth motion) does contribute to force generation in polymeric liquids (e.g. airway mucus.) 

In previous work, it was shown that viscoclastic stress systematically lead to a decrease of the swimming performance 



for swimmers who deform their flagella in a wave- like fashion 2J, |25| . [26[ . The calculations above display explicitly 
a class of body deformation where stresses arising from complex fluids act the other way, i. e. actually increase the 
forces generated. In the case of a general body deformation, both are in general possible and the beneficial vs. 
detrimental nature of viscoelastic stresses cannot be established a priori. It is also worth pointing out that if the 
flapper was flexible, it would deform as a result of fluid forces, and the resulting non-reciprocal shapes would lead to 



force generation even in the Newtonian limit [531 . |54 1 . 



Let us now estimate the order of magnitude of the viscoelastic forces reported here, and consider for simplicity the 
case of a single flapper near a wall (Fig. [1^). In that case, the result of Eqs. ((32|) leads to the order of magnitude 
estimate for the time-averaged force 



2 Aw rjs 



As an example of biological polymeric fluid, we consider cervical mucus |ii|,[i3,[iiQ,[i3, |16| . Since the viscosity of 

cervical mucus is at least two orders of magnitude that of water, we see that a large-amplitude flapping motion (e ~ 1) 
at frequency u ^ 1/X leads to |(^f2}| ^ ij^- Remarkably, stresses on the order of those that would be generated in a 
purely viscous fluid with viscosity 77 and sheared at a rate lo can be generated on average inside the polymeric fluid. 
For cervical mucus, we have rj ~ 10^^ — 10 Pas. Since A ^ 1 — 100 we would get uj ^ 10^^ — 1 s~^, leading to stresses 
on the order |(i5f2)| ^ 10~^ — 10 Pa, similar to the range of stresses applied for example by a ciliary array on an 
overlying water-like fluid. The effect reported in our paper could be therefore significant in a biological setting. 

A second important feature is to estimate how easy it would be to measure these viscoelastic forces in an experiment. 
Indeed, the forces on the flapper also include an unsteady component, which is of order e, and (although it averages 
to zero) might dominate the time-averaged force during experimental measurements. The amplitude of this unsteady 
term, in dimensional form, can be inferred from Eqs. (|18p - PD|) . and we see that 

When Xlu ^ 1, and e ~ 1 we obtain the estimate \6{i \ rjuj. Therefore for large-amplitude motion and a flapping rate 
comparable to the inverse relaxation time of the fluid, the unsteady forces on the flapper are expected to be of the 
same order as their time- average \S{i \ ^ |(<^f2)| ^ V-^- In that case, the breakdown of the scallop theorem would be 
significant, and could realistically be measured in an experiment. 

Finally, let us address the implication of our results to the locomotion of microorganisms. Our calculations assume 
the swimmer to be tethered, and therefore the analysis would have to be carried out separately in the case of free- 
swimming. This would also require to consider a finite-size body. Qualitatively, the time-averaged forces obtained 
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above would lead to propulsion. The organism would swim at a speed found by balancing the propulsive forces, 
Eq. (|47|) . with a viscoelastic drag. In theory, given that swimmers with arbitrarily complex geometries could be 
devised, locomotion could occur with arbitrarily complex trajectories. A systematic study of viscoelastic swimming 
for finite-size bodies will be reported in future work. 



C. Generalization to other fluids. 

Our calculations were performed for an Oldroyd-B fluid, which, in steady-shear, has a constant viscosity and first 
normal stress coefficient, and no second-normal stress differences. In previous work on the swimming of a sheet in a 
viscoelastic fluid 2J] , the Oldroyd-B results were extended to a more general class of Oldroyd-like models ( Johnson- 
Segalman-Oldroyd), to the (nonlinear) Giesekus model, both of which display shear-thinning behavior for the viscosity 
and both normal stress coefficients in steady shear. A more realistic polymeric model was also considered, FENE-P, 
which remains physically accurate for large values of the Deborah numbers. For the FENE-P model, the viscosity 
and first normal stress coefficient are shear-thinning, and the second normal stress coefficient is zero. 

Following the calculations in the appendices of Ref . [31 , it is straightforward to extend the flapping calculations 
above to these models. Our results, Eqs. ((32|) . (l33|) and ((46|), turn out to be unchanged (not reproduced here), and 
therefore appear to be quite general. 



D. Force vs. flow generation. 

An interesting feature of our asymptotic results is that the net forces occur at order O(e^). However, there is no 
time-averaged flow of the polymeric fluid at this order. The same result is obtained for all the constitutive models 
considered. A detailed (but lengthy) analysis reveals that a time-averaged flow actually occurs at order 0{e*), and 
will be reported in future work. How can we reconcile the fact that there could be a net force without a flow? In 
fact, a similar situation arises in the example of the cone-and-plate rheometer proposed above. Under the oscillatory 
motion of the rheometer, not net flow is generated, but time-averaged forces are arising. 

E. Asymptotic details. 

I — I 

As in the classic work by HK Moffatt 40], the calculations here should be thought of as similarity solution, valid 
away from the microscopic cuttoff length, tq at which the pressure divergence is regularized, but close enough to the 
hinge point that inertial effects can be neglected, i.e. eusr'^/h' ^ 1 or r ^ y/iT/euj, as well as edge effects for finite-size 
fiappers. 



15 



Acknowledgments 

Wc thank Henry Fu and Tom Powers for useful discussions. This work was supported in part by the National 
Science Foundation (Grants No. CTS-0624830 and CBET-0746285). 



[1] J. Gray. Animal locomotion. Norton, London, 1968. 

[2] R. M. Alexander. Principles of Animal locomotion. Princeton University Press, Princeton, NJ, 2003. 
[3] J. Lighthill. Mathematical Biofluiddynamics. SIAM, Philadelphia, 1975. 

[4] J. Lighthill. Flagellar hydrodynamics - The John von Neumann lecture, 1975. SIAM Rev., 18:161-230, 1976. 

[5] C. Brennen and H. Winet. Fluid mechanics of propulsion by cilia and flagella. Ann. Rev. Fluid Mech., 9:339-398, 1977. 

[6] S. Childress. Mechanics of Swimming and Flying. Cambridge Universtity Press, Cambridge U.K., 1981. 

[7] D. Bray. Cell Movements. Garland Publishing, New York, NY, 2000. 

[8] L. J. Fauci and R. Dillon. Biofluidmechanics of reproduction. Ann. Rev. Flmd Mech., 38:371-394, 2006. 
[9] S. S. Suarez and A. A. Pacey. Sperm transport in the female reproductive tract. Hunnan Reprod. Update, 12:23-37, 2006. 
[10] A. I. Yudin, F. W. Hanson, and D. F. Katz. Human cervical mucus and its interaction with sperm - A fine structural view. 

Biol. Reprod., 40:661-671, 1989. 
[11] D. P. Wolf, L. Blasco, M. A. Khan, and M. Litt. Human cervical mucus. I. Rheologic characteristics. Fertility Sterility, 
28:41-46, 1977. 

[12] D. P. Wolf, L. Blasco, M. A. Khan, and M. Litt. Human cervical mucus. II. Changes in viscoelasticity during ovulatory 

menstrual cycle. Fertility Sterility, 28:47-52, 1977. 
[13] D. P. Wolf, J. Sokoloski, M. A. Khan, and M. Litt. Human cervical mucus. HI. Isolation and characterization of Theologically 

active mucin. Fertility Sterility, 28:53-58, 1977. 
[14] D. P. Wolf, L. Blasco, M. A. Khan, and M. Litt. Human cervical mucus. IV. Viscoelasticity and sperm penetrability during 

ovulatory menstrual-cycle. Fertility Sterility, 30:163-169, 1978. 
[15] D. P. Wolf, J. E. Sokoloski, and M. Litt. Composition and function of human cervical mucus. Biochimica Biophysica Acta, 

630:545-558, 1980. 

[16] P. Y. Tam, D. F. Katz, and S. A. Berger. Nonlinear viscoelastic properties of cervical mucus. Biorheol., 17:465-478, 1980. 
[17] J. B. Shukla, B. R. P. Rao, and R. S. Parihar. Swimming of spermatozoa in cervix - Effects of dynamical interaction and 

peripheral layer viscosity. J. Biomech., 11:15-19, 1978. 
[18] D. F. Katz, R. N. Mills, and T. R. Pritchett. Movement of human spermatozoa in cervical mucus. J. Reprod. Fert., 

53:259-265, 1978. 

[19] D. F. Katz and S. A. Borgor. Flagellar propulsion of human sperm in cervical mucus. Biorheol., 17:169-175, 1980. 
[20] D. F. Katz, T. D. Bloom, and R. H. Bondurant. Movement of bull spermatozoa in cervical mucus. Biol. Reprod., 25:931-937, 
1981. 

[21] R. Rikmenspoel. Movements and active moments of bull sperm flagella as a function of temperature and viscosity. J. Exp. 
Biol., 108:205-230, 1984. 

[22] S. Ishijima, S. oshio, and H. mohri. Flagellar movement of human-spermatozoa. Gam. Res., 13:185-197, 1986. 

[23] S. S. Suarez and X. B. Dai. Hyperactivation enhances mouse sperm capacity for penetrating viscoelastic media. Biol. 

Reprod., 46:686-691, 1992. 
[24] E. Lauga. Propulsion in a viscoelastic fluid. Phys. Fluids, 19:083104, 2007. 

[25] H. Fu, T. R. Powers, and C. W. Wolgemuth. Theory of swimming filaments in viscoelastic media. Phys. Rev. Lett, 
99:258101, 2008. 

[26] H. Fu, C. W. Wolgemuth, and T. R. Powers. Swimming speeds of filaments in nonlinearly viscoelastic fluids, preprint, 



16 



2008. 

[27] R. B. Bird, R. C. Armstrong, and O. Hassager. Dynamics of Polymeric Liquids. Second Edition. Vol. 1: Fluid Mechanics. 

Wiley-Interscience, New York, NY, 1987. 
[28] R. B. Bird, C. F. Curtiss, R. C. Armstrong, and O. Hassager. Dynamics of Polymeric Liquids. Second Edition. Vol. 2: 

Kinetic Theory. Wiley-Interscience, New York, NY, 1987. 
[29] M. Doi and S. F. Edwards. The Theory of Polymer Dynamics. Oxford University Press, Oxford, U.K., 1988. 
[30] R. G. Larson. The Structure and Rheology of Complex Fluids. Oxford Universtity Press, Oxford, U.K., 1999. 
[31] E. M. Purcell. Life at low Reynolds number. Am. J. Phys., 45:3-11, 1977. 

[32] S. Childress and R. Dudley. Transition from ciliary to flapping mode in a swimming mollusc: Flapping flight as a bifurcation 

in Rcu,. J. Fluid Mech., 498:257-288, 2004. 
[33] N. Vandenberghe, J. Zhang, and S. Childress. Symmetry breaking leads to forward flapping flight. J. Fluid Mech., 

506:147-155, 2004. 

[34] S. Albon and M. ShcUcy. Coherent locomotion as an attracting state for a free flapping body. Proc. Natl. Acad. Sci. USA, 
102:11163-11166, 2005. 

[35] E. Lauga. Continuous breakdown of Purcell's scallop theorem with inertia. Phys. Fluids, 19:061703, 2007. 
[36] M. A. Sleigh, J. R. Blake, and N. Liron. The propulsion of mucus by cilia. Am. Rev. Resp. Dis., 137:726-741, 1988. 
[37] G. I. Taylor. Similarity solutions of hydrodynamic problems, pages 21-28. Aeronautics and Astronautics. Pergamon Press, 
Oxford, U. K., 1960. 

[38] G. K. Batchelor. An Introduction to Fluid Dynamics. Cambridge University Press, Cambridge, UK, 1967. 
[39] C. P. Hills and H. K. Moffatt. Rotary honing: a variant of the taylor paint-scraper problem. J. Fluid Mech., 418:119-135, 
2000. 

[40] H. K. Mofi^att. Viscous and resistive eddies near a sharp corner. J. Fluid Mech., 18:1-18, 1964. 

[41] H.K. Moffatt and B.R. Duffy. Local similarity solutions and their limitations. J. Fluid Mech., 96:299-313, 1980. 

[42] C. Chang and W. Schowalter. Flow near an oscillating cylinder in dilute viscoelastic fluid. Nature, 252:686-688, 1974. 

[43] P. W. James. Elastico-viscous flow around a circular-cylinder executing small amplitude, high-frequency oscillations. J. 

Non-Newt. Fluid Mech., 2:99-107, 1977. 
[44] S. Rosenblat. Reversal of secondary flow in non-newtonian fluids. J. Fluid Mech., 85:387-400, 1978. 
[45] G. Bohme. On steady streaming in viscoelastic liquids. J. Non-Newt. Fluid Mech., 44:149-170, 1992. 

[46] K. C. Bagchi. Flow of non-Newtonian fluid contained in a fixed circular cylinder due to longitudinal oscillation of a 

concentric circular cylinder. Appl. Sci. Res., 16:131-140, 1966. 
[47] K. R. Frater. Acoustic streaming in an elastico-viscous fluid. J. Fluid Mech., 30:689-697, 1967. 

[48] K. R. Frater. Drag on a circular cylinder oscillating in an elastico-viscous fluid. Z. Angew. Math. Mech., 19:510-512, 1968. 
[49] C. Goldstein and W. R. Schowalter. Further studies of fluid nonlinearity - Orthogonal rheometer and oscillating sphere. 

Trans. Soc. RheoL, 19:1 19, 1975. 
[50] C. F. Chang. Boundary-layer analysis of oscillating cylinder flows in a viscoelastic liquid. Z. Angew. Math. Mech., 28:283- 
288, 1977. 

[51] C. F. Chang and W. R. Schowalter. Secondary flow in the neighborhood of a cylinder oscillating in a viscoelastic fluid. J. 
Non-Newt. Fluid Mech., 6:47-67, 1979. 

[52] M. L. Roper. Symmetry Breaking and Un-breaking in Microhydrodynamical Systems: Swimming, Pumping and Bio- 
ballistics. PhD thesis. School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 2007. 
[53] K. E. Machin. Wave propagation along flagella. J. Exp. Biol, 35:796-806, 1958. 

[54] C. H. Wiggins and R. E. Goldstein. Flexive and propulsive dynamics of elastica at low Reynolds number. Phys. Rev. Lett., 
80:3879 - 3882, 1998. 



